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Abstract 


Background: Following the initial isolation of porcine deltacoronavirus (PDCoV) from pigs with diarrheal disease in 
the United States in 2014, the virus has been detected on swine farms in some provinces of China. To date, little is 
known about the molecular epidemiology of PDCoV in southern China where major swine production is operated. 


Results: To investigate the prevalence of PDCoV in this region and compare i 
caused by porcine epidemic diarrhea virus (PEDV), transmissible gastroenteritis 
group C (Rota ©), 390 fecal samples were collected from swine of various ages 
Fecal samples were tested by reverse transcription-PCR (RT-PCR) that targeted 
PDCoV was detected exclusively from nursing piglets with an overall prevalence of approxima 


suckling and fattening piglets. Interestingly, all of PDCoV-posi 
a low detection rate, PDCoV emerged in each province/region of southern Ch 


5/390) or Rota C (1.28 %, 6/390), there were highly detection 
five PDCoV-positive piglets were co-infected by PEDV. Furth 
(N) gene sequences of PDCoVs revealed that currently circ 
to other Chinese strains of PDCoVs than to those reported 


ts activity to other enteric disease of swine 
coronavirus (TGEV), and porcine rotavirus 
from 15 swine farms with reported diarrhea. 
PDCoV, PEDV, TGEV, and Rota C, respectively. 
te 1.28 % (5/390), not in 

tive samples were from 2015 rather than 2012-2014. Despite 
ina. In addition, compared to TGEV (1.54 %, 
rates of PEDV (22.6 %, 88/390) in those samples. Notably, all 
ermore, phylogenetic analysis of spike (S) and nucleocapsid 


ulating PDCoVs in southern China were more closely related 


in United States, South Korea and Thailand. 


Conclusions: This study demonstrated that PDCoV was present in southern China despite the low prevalence, and 


supported an evolutionary theory of geographical clustering of PDCoVs. 
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Background 

Before 2012, the subfamily Coronavirinae included three 
genera (Alphacoronavirus, Betacoronavirus and Gamma- 
coronavirus). However, in 2012, an emerging genus, Delta- 
coronavirus, was found in many animal species including 
swine from Hong Kong [1]. At present, more than five 
different coronaviruses have been described in swine pop- 
ulations. Among them, porcine epidemic diarrhea virus 
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(PEDV), transmissible gastroenteritis virus (TGEV), and 
porcine respiratory coronavirus (PRCV) belong to the genus 
Alphacoronavirus; meanwhile, porcine hemagglutinating en- 
cephalomyelitis virus (PHEV) and Porcine deltacoronavirus 
(PDCoV) are assigned to the genus Betacoronavirus and 
the genus Deltacoronavirus, respectively [1]. Numerous 
studies have shown that more than half of porcine corona- 
viruses (including PDCoV) were enteropathogenic and 
caused acute diarrhea and vomiting in pigs, which resulted 
in huge economic losses for the global swine industry [2-6]. 

Currently, PDCoV has been reported in Hong Kong, 
North America, Mexico, South Korea, Thailand and some 
provinces of China [1, 7-19]. Despite recent progress, 
little is known about the prevalence and epidemiology of 
PDCoV in southern China (including Guangdong 
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province, Hainan province, and Guangxi autonomous 
region), where major swine production is operated. 
Therefore, the aim of this study was to investigate the 
prevalence and sequence properties of PDCoV in this 
region. 


Methods 

Sampling 

A total of 390 fecal samples (Table 1) were collected 
from 15 commercial swine farms with reported diarrhea 
in southern China. Farms A, D, E, J-O were from 
Guangdong province, Farms B, F, G were from Hainan 
province, and farms C, H, I were the Guangxi autono- 
mous region. Farms A-I derived 30 samples with the 
following arrangement: ten samples from suckling 
piglets (<3 weeks old), ten samples from nursing piglets 
(between 3 and 9 weeks old), and ten samples from 
fattening piglets (>9 weeks old). The samples of farms 
A-I were collected between July and August 2015 and 
stored at —-80 °C until further use. However, the samples 
of farms J-O were archived samples from 2012 to 2014. 
Prior to viral RNA extraction, fecal samples were diluted 
one time using Phosphate Buffered Saline (PBS) (pH: 
7.4). The supernatants were then collected by centrifuga- 
tion at 5000 x g for 5 min. 200 ul of clarified superna- 
tants was used to extract viral RNA following the 
manufacturer’s recommendations (Axygen Scientific 
Inc.). RNA samples were stored at -80 °C until further 
analysis. 


Reverse transcription polymerase chain reaction (RT-PCR) 
detection 

To detect PDCoV genome in collected fecal samples, 
the previously reported RT-PCR primers (41 F: 5’-TTT 
CAGGTGCTCAAAGCTCA-3’ and 735R: 5’-GCGAAA 
AGCATTTCCTGAAC-3’) targeting the nucleocapsid 
(N) gene with reaction conditions (50 °C for 30 min 
and 95 °C for 15 min for the reverse transcription reac- 
tion, followed by 40 cycles of PCR amplification at 95 °C 
for 15 s, 55 °C for 45 s, and 72 °C for 1 min, with a final 
extension at 72 °C for 7 min) were used [15]. In addition, 
molecular detection of the three diarrhea-related enteric 
viruses (Porcine epidemic diarrhea virus, PEDV; Porcine 
transmissible gastroenteritis virus, TGEV; Porcine rota- 
virus group C, Rota C) was performed in accordance with 
previous methods [20-22] for further evaluation of the 
possible co-infection status with PDCoV in investigated 
pig samples. 


Amplification of the spike (S) and N genes 

To perform in-depth sequence comparison and phylogen- 
etic analysis with known reference sequences (Additional 
file 1: Table S1), the complete spike (S) and N genes of 
PDCoV-positive samples were amplified according to 
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previously published methods [8]. For amplification of the 
full-length S and N genes, previously reported RT-PCR 
primers (PDCoV-SF2: 5’°-AGCGTTGACACCAACCTA 
TT-3’ and PDCoV-SR2: 5’-TCGTCGACTACCATTCCT 
TAAAC-3’; PDCoV-NF1 : 5’°-CCATC GCTCCAAG TC 
ATTCT-3’ and PDCoV-NR1: 5-TGGGTGGGTTTAA 
CAGACATAG-3’) were used. PCR was carried out at 
50 °C 30 min and 95 °C for 5 min, followed by 40 cycles 
of 98 °C for 10 s, 55 °C for 15 s, and 68 °C for 5 and 2 min 
for S and N genes, respectively; final extension was per- 
formed at 68 °C for 15 min. Positive amplicons were 
cloned into the pGM-19 T vector (Tiangen Inc. Beijing). 
Furthermore, all positive recombinant plasmids were sub- 
mitted to a sequencing company (The Beijing Genomics 
Institute, BGI) and sequenced at least three times. Five S 
gene sequences and five N gene sequences were obtained 
(Additional file 1: Table S1), and have been submitted 
to GenBank database (accession numbers KU204694- 
KU204701, KX534090-KX534091). 


Phylogenetic analysis of the S and N genes 

Sequence alignment analysis was performed using the 
Clustal W program implemented in DNAStar software. A 
phylogenetic tree was then constructed by the neighbor- 
joining method using the Molecular Evolutionary Genetics 
Analysis (MEGA) software version 5.1 with 1000 bootstrap 
replications set at 1000. Moreover, the possible recombin- 
ation event was evaluated in the S and N genes by recom- 
bination detection program (RDP) 3.34 software. 


Results 

PDCoV detection 

A total of 390 pig fecal samples, collected from 15 swine 
farms with reported diarrhea in southern China, were 
assessed for the presence of PDCoV and other viral en- 
teric pathogens (PEDV, TGEV, and Rota C) by RT-PCR. 
As summarized in Table 1, the PDCoV genome was de- 
tected in specimens from 4 of 15 swine farms. Interest- 
ingly, the PDCoV genome was detected only in nursing 
piglets, and was absent in suckling piglets and fattening 
pigs. Although PDCoV was detected in each province/ 
region of southern China, its overall prevalence in the 
investigated pigs of various age groups (m= 390) was 
relatively low (5/390, 1.28 %). The positive rate could be 
higher if only nursing piglets were included (5/150, 
3.33 %). In contrast, the prevalence of PEDV, another por- 
cine coronavirus causing epidemic diarrhea, was relatively 
higher (22.6 %, 88/390). In addition, PEDV was different 
from PDCoV in that it distributed similarly between nurs- 
ing (36/150, 24 %) and suckling piglets (47/150, 31.3 %). 
Five fattening pigs from farms A, D, F and I were also 
tested positive for the PEDV genome. We also examined 
whether pigs with diarrhea harbored other enteric viruses 
such as TGEV and Rota C. Our results showed that the 
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Table 1 Sample information and RT-PC 
farms in southern China 


R detection results of four diarrhea-associated viruses in pigs of various ages from 15 swine 
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Farm Pigs (Age) Number (n) PDCoV PEDV TGEV Rota C 
A Suckling piglets (<3 weeks old) 0 0/10 2/10 0/10 0/10 
ursery pigs (>3 weeks old, <9 weeks old) 0) 2/10 2/10 0/10 0/10 
Fattening pig (>9 weeks old 0 0/10 1/10 0/10 0/10 
B Suckling piglets (<3 weeks old) 0 0/10 3/10 0/10 0/10 
ursery pigs (>3 weeks old, <9 weeks old) 0) 1/10 2/10 0/10 0/10 
Fattening pig (>9 weeks old 0 0/10 0/10 0/10 0/10 
C Suckling piglets (<3 weeks old) 0) 0/10 3/10 0/10 0/10 
ursery pigs (>3 weeks old, <9 weeks old) 0) 1/10 3/10 0/10 0/10 
Fattening pig (>9 weeks old 0) 0/10 0/10 0/10 0/10 
D Suckling piglets (<3 weeks old) 0) 0/10 2/10 0/10 0/10 
ursery pigs (>3 weeks old, <9 weeks old) 0) 0/10 1/10 0/10 2/10 
Fattening pig (>9 weeks old 0 0/10 1/10 1/10 0/10 
E Suckling piglets (<3 weeks old) 0) 0/10 1/10 0/10 0/10 
ursery pigs (>3 weeks old, <9 weeks old) 0 1/10 2/10 0/10 0/10 
Fattening pig (>9 weeks old 0) 0/10 0/10 0/10 0/10 
F Suckling piglets (<3 weeks old) 0) 0/10 3/10 0/10 0/10 
ursery pigs (>3 weeks old, <9 weeks old) 0 0/10 4/10 0/10 0/10 
Fattening pig (>9 weeks old 0 0/10 1/10 2/10 0/10 
G Suckling piglets (<3 weeks old) 0 0/10 5/10 0/10 0/10 
ursery pigs (>3 weeks old, <9 weeks old) 0) 0/10 3/10 0/10 0/10 
Fattening pig (>9 weeks old 0) 0/10 0/10 0/10 0/10 
H Suckling piglets (<3 weeks old) 0 0/10 7/10 0/10 0/10 
ursery pigs (>3 weeks old, <9 weeks old) 0) 0/10 2/10 0/10 0/10 
Fattening pig (>9 weeks old 0) 0/10 0/10 0/10 0/10 
Suckling piglets (<3 weeks old) 0) 0/10 4/10 0/10 0/10 
ursery pigs (>3 weeks old, <9 weeks old) 0) 0/10 2/10 0/10 0/10 
Fattening pig (>9 weeks old 0 0/10 2/10 0/10 0/10 
J Suckling piglets (<3 weeks old 0 0/10 4/10 0/10 0/10 
ursery pigs (>3 weeks old, <9 weeks old) 0 0/10 3/10 0/10 0/10 
Suckling piglets (<3 weeks old 0 0/10 5/10 0/10 2/10 
ursery pigs (>3 weeks old, <9 weeks old) 0 0/10 5/10 2/10 0/10 
L Suckling piglets (<3 weeks old 0 0/10 2/10 0/10 0/10 
ursery pigs (>3 weeks old, <9 weeks old) 0 0/10 0/10 0/10 0/10 
Suckling piglets (<3 weeks old 0 0/10 2/10 0/10 0/10 
ursery pigs (>3 weeks old, <9 weeks old) 0 0/10 3/10 0/10 0/10 
Suckling piglets (<3 weeks old 0 0/10 2/10 0/10 1/10 
ursery pigs (>3 weeks old, <9 weeks old) 0 0/10 1/10 0/10 0/10 
O Suckling piglets (<3 weeks old 0 0/10 2/10 0/10 0/10 
ursery pigs (>3 weeks old, <9 weeks old) 0 0/10 3/10 1/10 0/10 

Total 390 5/390 88/390 6/390 5/390 


low detection rates (1.54 %, 5/390 for TGEV vs 1.28 %, 
6/390 for Rota C) of the two pathogens were present in 
those pig samples. Intriguingly, co-infection of pigs by 


PDCoV and PEDV was observed (Table 1). All PDCoV 
positive nursing piglets were also tested positive for PEDV, 
thereby indicating a 100 % co-infection rate. 
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Sequence comparison and phylogenetic analysis of the S 
gene of PDCoVs 

The full-length sequences of S genes in five PDCoV- 
positive samples from the four different farms were 
amplified and designated provisionally CH/GD01/2015, 
CH/GD02/2015, CH/GD03/2015, CH/HN01/2015, and 
CH/GX01/2015, respectively. Sequencing results showed 
that they were composed of 3480 nucleotides (nt). Com- 
pared to all published American and individual Asian 
strains (including HKU15-44, CHN-AH-2004, KNU14- 
04, PDCoV/Swine/Thailand/S5015L/2015 and PDCoV/ 
Swine/Thailand/S5011/2015), a 3-nt deletion in the S 
gene was identified in the five current PDCoV strains, six 
other Chinese viral strains and one Thai viral strain reported 
previously (Additional file 1: Table $1) [5, 17, 18]. Sequence 
alignment results revealed that CH/GD01/2015, CH/GD02/ 
2015, CH/GD03/2015, CH/HN01/2015 and CH/GX01/ 
2015 had >99.9 % homology in the S gene nucleotide 
sequence, indicating these five viral strains evolved from the 
same ancestor (Table 2). Further analysis showed that the 
five viruses reported in this study had the highest nucleotide 
identity (98.9 to 99.5 %) with a Jiangxi strain (PDCoV/ 
CHJXNI2/2015) isolated from Jiangxi province bordering 
with the northern region of Guangdong province, and 
possessed the lowest nucleotide similarity (95.4 to 95.9 %) 
with PDCoV/Swine/Thailand/S5011/2015 and PDCoV/ 
Swine/Thailand/S5015L/2015, isolated from Thailand 
(Table 2). From the above data, phylogenetic analysis of the 
S gene showed that the current PDCoVs circulating in 
southern China were most closely related to other Chinese 
PDCovV isolates than to those isolated previously from 
USA, South Korea and Thailand (Fig. 1). In addition, in the 
S gene, any possible recombinant events were not detected 
among those PDCoV strains. 


Sequence comparison and phylogenetic analysis of the N 
gene of PDCoVs 

Similarly, N gene sequences of all five PDCoV-positive sam- 
ples were identified as 1029 nt in size. Sequence alignment 
results suggested that there was no deletion or insertion in 
N gene regions (Additional file 1: Table $1). Consistent with 
the data of S gene, multiple sequence alignment results of 
N gene showed that CH/GD01/2015, CH/GD02/2015, CH/ 
GD03/2015, CH/HNO01/2015 and CH/GX01/2015 had the 
highest nucleotide homology (99.1 to 99.7 %) with PDCoV/ 
CHJXNI2/2015, and the lowest nucleotide homology (96.9 
to 97.1 %) with PDCoV/Swine/Thailand/S5011/2015 
and PDCoV/Swine/Thailand/ S5015L/2015 (Table 2). In 
addition, in the phylogenetic tree based on N gene, PDCoVs 
were divided into three main branches (Chinese branch, 
American branch and Thai branch). The five viral isolates 
reported from this work were clustered together within the 
Chinese branch (Fig. 2). Moreover, there were no any 
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possible recombinant events occurring at the N gene of 
those PDCoV strains. 


Discussion 

Epidemiology of PDCoVs 

PDCoV was first identified by Deltacoronavirus specific- 
PCR in rectal swabs of pigs (10.1 %, 17/169) with un- 
known healthy status in Hong Kong [1]. Then, PDCoV 
emerged in United States, China, South Korea and 
Thailand [5, 17, 18]. In most of studies, excluding PEDV, 
TGEV and porcine rotavirus, PDCoV, as an important 
enteric pathogen, was detected in clinical samples from 
pigs with diarrhea [9-12, 23]. In addition, it was con- 
firmed experimentally that less than two-week old piglets 
were susceptible to PDCoV, which caused mild to moder- 
ate diarrhea as well as macroscopic and microscopic lesions 
in small intestines of conventional piglets (5-day-old), and 
severe diarrhea, vomiting, fecal shedding of virus, and se- 
vere atrophic enteritis in gnotobiotic pigs (11- to 14-day- 
old) [2, 3]. The data further confirmed that PDCoV were 
enteropathogenic in pigs. Meanwhile, PEDV or rotavirus 
showed higher detection rates in PDCoV-positive samples 
compared with other TGEV and rotavirus [8, 13-15, 24]. 
As shown above, co-infection of PDCoV and PEDV oc- 
curred in nursing piglets (Table 1), indicating that the 
diarrhea-related pathogens were quite complex clinically 
and not easy to control in the field. Moreover, in the two 
recent studies, PDCoV was shown to have higher infectiv- 
ity in sows with diarrhea (81.0 %, 34/42) than nonclinical 
counterparts (23.5 %, 4/17) [8, 15], which might imply that 
sows often carry PDCoV. And further, it could result in the 
transmission of PDCoV from sows to the foetus and even 
newborn piglets, although the pathogenesis mechanism of 
PDCoV remains unclear. 

To further understand the origin of PDCoV, some retro- 
spective studies were made using PCR and enzyme-linked 
immunosorbent assay (ELISA) [24—26]. In Dong et al. [24] 
study, 2 of 6 samples collected from Anhui Province of 
China in 2004 were positive for PDCoV, up to now, it was 
the most ancient time for the detection of PDCoV in 
China. Meanwhile, PDCoV could date back to August 
2013 in United States, where only 3 PDCoV-samples were 
detected using PCR in archived samples [25]. As for 
serology of PDCoV, anti-PDCoV IgG antibodies could 
date back to 2010 using an indirect anti-PDCoV IgG 
ELISA based on the putative S1 portion of the spike pro- 
tein [26]. The above studies indicated that PDCoV could 
have circulated in China at least since 2004 and in United 
States since 2010. Maybe, due to limted samples in the 
present study, we did not detect PDCoV in pig samples col- 
lected from Guangdong province between 2012 and 2014. 
Although Asian leopard cat coronavirus (GenBank 
accession no. EF584908) was close to PDCoV in the phylo- 
genetic trees (Figs. 1 and 2), in the future, more 


Table 2 Nucleotide similarities (%) of S and N genes of our five PDCoV isolates and other reported PDCoVs and coronaviruses 


S gene N gene 
GenBank No. KU204694 KU204695 KU204696 KU204697 KX534090 KU204698 KU204699 KU204700 KU204701 KX534091 
Strain CH/GD01/ CH/GD02/ CH/HNO1/ CH/GX01/ CH/GD03/ CH/GD01/ CH/GD02/ CH/HNO1/ CH/GX01/ CH/GD03/ 
2015 2015 2015 2015 2015 2015 2015 2015 2015 2015 
KR131621 PDCoV/CHJXNI2/2015 99.5 98.9 98.9 99.1 99.2 99.7 99.1 99.1 993 99.4 
KP757892 CHN-JS-2014 99.1 98.5 98.5 98.7 98.7 99,3 99,3 99.3 99.5 99.0 
KP757891 CHN-HB-2014 98.9 98.4 98.4 98.7 98.5 98.8 98.6 98.6 98.8 98.5 
JQ065042 HKU15-44 98.6 98.2 98.2 98.5 98.3 98.7 98.5 98.5 98.7 98.4 
JQ065043 HKU15-155 98.5 98.2 98.2 98.4 98.2 98.7 98.5 98.5 98,7 98.4 
KT021234 CH/SXD1/2015 98.4 979 97.9 98.1 98.0 99.0 98.8 98.8 99.0 98.7 
KT266822 CH/Sichuan/S27/2012 98.4 98.0 98.0 98.2 98.0 99.2 99.0 99.0 99.2 98.9 
KM820765 U14-04 98.4 97.9 97.9 98.2 98.1 98.8 98.6 98.6 98.8 98.5 
KJ620016 16148 98.4 97.9 97.9 98.1 98.0 98.8 98.6 98.6 98.8 98.5 
KJ584360 N3092 -a -a -a -a -a 98.8 98.6 98.6 98.8 98.5 
KJ584358 PA3148 98.4 98.0 98.0 98.2 98.0 98.8 98.6 98.6 98.8 98.5 
KJ584357 4813 98.4 97.9 97.9 98. 98.0 98.8 98.6 98.6 98.8 98.5 
KJ584355 L2768 98.4 97.9 97.9 98. 98.0 98.8 98.6 98.6 98.8 98.5 
KT381613 OH11846 98.4 97.9 97.9 98. 98.0 98.8 98.6 98.6 98.8 98.5 
KJ601779 PDCOV/USA/IIlinois 98.4 979 97.9 98. 98.0 98.5 98.3 98.3 98.5 98.3 
36/2014 
KJ481931 PDCOV/USAMIlinois 98.4 979 97.9 98. 98.0 98.5 98.3 98.3 98.5 98.3 
21/2014 
KJ769231 OhioCVM1/2014 98.4 978 97.8 98. 98.0 98.3 98.2 98.2 98.3 98. 
KJ601777 PDCOV/USA/IIlinois 98.4 978 978 98. 98.0 98.7 98.5 98.5 98.7 98.4 
33/2014 
KJ584359 E3579 98.4 978 978 98. 98.0 98.7 98.5 98.5 98.7 98.4 
KJ584356 SD3424 98.4 97.8 978 98. 98.0 98.5 98.3 98.4 98.5 98.3 
KJ462462 OH1987 98.4 978 978 98. 98.0 98.7 98.5 98.5 98,7 984 
KJ601778 PDCOV/USATIIlinois 98.3 978 97.8 98.0 98.0 98.7 98.5 98.5 98,7 984 
134/2014 
KP995358 OH-FD22 98.3 978 978 98.0 98.0 -b -b -b -b -b 
KJ601780 PDCoV/USA/Ohio 98.3 978 97.8 98.0 98.0 98.7 98.5 98.5 98,7 984 
137/2014 
KJ569769 IN2847 98.3 978 978 98.1 98.0 98.8 98.6 98.6 98.8 98.5 
KJ567050 8734/USA-IA/2014 98.3 979 97.9 98.1 98.0 98.7 98.5 98.5 98,7 984 
KP981395 98.3 97.8 978 98.0 97.9 98.7 98.5 98.5 98.7 98.4 
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Table 2 Nucleotide similarities (%) of S and N genes of our five PDCoV isolates and other reported PDCoVs and coronaviruses (Continued) 


USA/IL/2014/026 
PDV_P11 
012168 ichigan/8977/2014 98.3 978 978 98.0 97.9 98.7 98.5 98.5 98.7 98.4 

P757890 CHN-AH-2004 98.0 97.5 9755 97.8 97.7 98.1 97.9 97.9 98.1 97.8 
U051641 PDCoV/Swine/Thailand/ 95.8 95.6 95.6 95.9 95.4 97.1 96.9 96.9 97.1 96.8 

$5011/2015 
U051649 PDCoV/Swine/Thailand/ 95.8 95.5 95.5 95.9 95.4 97.1 96.9 96.9 97.1 96.8 

$5015L/2015 
U984334 P23_15_TT_1115 95.9 95.8 95.6 95.9 95.5 97.8 97.8 97.8 98.0 97.5 
EF584908 Guangxi/F230/2006 97.5 97.1 97.1 97.3 97.0 98.0 978 97.8 98.0 97.7 
JQ065045 HKU17-6124 40.9 40.8 40.8 40.7 40.7 92.9 92.6 92.8 92.9 92.6 
FJ376621 HKU 12-600 40.7 41.1 41.1 41.0 40.6 76.4 75:7 76.2 76.3 76.3 
JQ065044 HKU16-6847 57.0 56.9 56.9 57.0 57.2 75.1 75.3 75.2 754 75.1 
FJ376619 HKU1 1-934 65.9 68.9 68.9 66.1 65.8 744 74.3 745 744 74.2 
FJ376620 HKU1 1-796 65.3 65.6 65.7 65.6 65.1 744 746 745 744 74.2 
KJ408801 OH1414 (PEDV) 38.2 38.3 38.3 38.1 378 3.6 3.2 3.5 3.6 3.6 
FJ755618 H16(TGEV) 38.1 38.1 38.1 38.0 37.9 28 49 75 5.0 28 
DQ811787 SU-1(PRCV) 374 25.7 25.7 25.7 37.3 73 7A Lo 74 7A 
BCU00735 ebus(bovine) 124 12.5 12.5 124 71 Zi 2. 2.1 2.1 2. 
AY654624 TJF(SARS) 6.8 6.9 6.9 68 6.9 8.8 88 24 8.8 8. 
DQ011855 VW572(PHEV) 21,9 21.9 21.9 21.9 218 2.1 2, 2.1 2.1 2: 
JF893452 YN(CIBV) 22.0 22.1 22.1 22.3 22.0 125 12.5 125 12.6 124 
NC_010800 MG10 (Turkey) 24.1 23.5 23.5 23.3 24.0 14.0 14.0 14.0 14.1 14.1 
NC_010646 SW1 (whale) 15.0 13.7 13.7 13.8 148 77 ra 77 77 77 
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3s gene of MN3092 was not complete; °S gene of OH-FD22 was not available 
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Strains/GenBank No.s Year Genus 
400 >- Mebus (BCU00735) * . 
vf L WWs72 (DOO11855) 1972 | Betacoronavirus 
100 7 
s[ La eon a }Gammacoronavirus 
m TJF (AY654624) 2003 | Alphacoronavirus 
a ee : | Gammacoronavirus 
-— OH1414 (KJ408801) 2014 
tS ogg SU aatz87) : } Alphacoronavirus 
400 + H16 (FJ755618) 1973 
HKU12-600 (FJ376621) 2007 
L HKU17-6124 (JQ065045) 2007 
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Fig. 1 Phylogenetic analysis based on the S gene of PDCoVs and other coronaviruses. Note: Those PDCoV strains from China, America, South 
Korea and Thiland were labelled using “m’, “e’, “@" and “a”, respectively. Moreover, PDCoV strains in this study were labelled using left arrows. 
The collection time was not available for those coronavirus strains labelled using star symbols 
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Fig. 2 Phylogenetic analysis based on the N gene of PDCoVs and other coronaviruses. Note: Those PDCoV strains from China, America, South 
Korea and Thiland were labelled using “m’, “e’, “@" and “a”, respectively. Moreover, PDCoV strains in this study were labelled using left arrows. 
The collection time was not available for those coronavirus strains labelled using star symbols 
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epidemiological surveys should be warranted to uncover 
the origin of PDCoV. 

At the territory of China, Southern China mainly 
includes Guangdong province, Hainan province and the 
Guangxi autonomous region. Although molecular 
detection of PDCoV was performed in these regions 
[23, 24], little information was available on PDCoV 
prevalence. In a study by Chen et al. [23], an overall 
positive-PDCoV rate of 23.4 % (15/64) was obtained in 
all samples collected from Guangdong, Shanxi and 
Hubei provinces. However, more detailed data of 
PDCoV was not available in Guangdong province. 
Meanwhile, in the study from Dong et al. [24], only 
four archived samples from the Guangxi autonomous 
region were examined, but all negative for PDCoV. In 
this study, we demonstrated that PDCoV circulated and 
was co-infected by PEDV on those swine farms in 
Guangdong province, Hainan province and the Guangxi 
autonomous region, which further contributed to the 
epidemiology of PDCoV in these regions despite the 
relatively low prevalence. 


Genetic diversity of PDCoVs 

The first two reported full-length PDCoV genome se- 
quences (HKU15-44 and HKU15-155) were 25, 437 nt 
and 25, 432 nt in length, respectively [1], and they had 
99.1 % nucleotide similarity with each other. Moreover, 
further sequence alignment showed a 3-nt (TAA) inser- 
tion in the S gene and a 3-nt (TTA) insertion in the 3’ 
untranslated region (UTR) of HKU15-44 [1, 5]. During 
the past 3 years, PDCoV-associated swine enteric disease 
was paid great attention in the major pig producing 
countries, especilly United States and China. Up to May 
2016, more than 30 complete PDCoV genome sequences 
were published in the GenBank database. All were 
generated in China and United States except for one 
sequence from South Korea and three sequences from 
Thailand [8-19]. The Korean strain, KNU14-04, had 25, 
422 nt in length, with similar genome features (a 3-nt 
insertion in the S gene with 3, 483 nt and a 3-nt inser- 
tion in the 3’ UTR, respectively) to all American strains 
and the Chinese strain (HKU15-44) [9]. Comparing 
complete genomes of the remaining Chinese strains to the 
American, Thai and Korean counterparts, CHN-HB-2014, 
CHN-JS-2014, PDCoV/CHJXNI2/2015, CH/Sichuan/S27/ 
2012, CH/SXD1/2015 and P23_15_TT_1115 only had the 
3-nt (AAT) deletion in the S gene (3, 480 nt) [14, 23, 24], 
while CHN-AH-2004 only had the 3-nt (TAA) deletion 
in the 3’ UTR [5, 24]. In the present study, 3-nt inser- 
tion was not found in UTR for our five obtained 
PDCoV (data not shown). Moreover, two additional 
unique features including a 6-nt (TTTGAA) deletion in 
the nonstructural protein (nsp) 2 region and a 9-nt 
(GCCGGTTGG) deletion in the nsp 3 region were also 
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found in CH/Sichuan/S27/2012 [16]. However, for Thai 
viral isolates, they owned one additional unique nucleotide 
(C) insertion in the 3’ UTR [17, 18]. The biological signifi- 
cance of these naturally occurring deletions or insertions 
in PDCoV biology and pathogenesis warrants further 
investigations. 

In this study, five S and five N gene sequences, respect- 
ively, were obtained to evaluate wherever genetic diversity 
of PDCoVs existed in southern China. Our results showed 
that these five S and five N gene sequences were more 
closely related to Chinese strains, and all clustered together 
in the phylogenetic tree (Table 2, Figs. 1 and 2). However, 
CH/GD01/2015 and CH/GD02/2015, reported in this 
study, originated from the same pig farm in Guangdong 
province, but had 48 nt and 12 nt differences in the S and 
N genes, respectively. The observed 48 nucleotide changes 
in the S gene made these viruses differ by 25 amino acid 
residues (Additional file 2: Table S2). For the N gene, the 
12 nucleotide changes among these viruses resulted in 3 
amino acid substitutions (Additional file 2: Table S2). 
Among them, 18 of 25 amino acid differences occurred at 
the first two-third parts of S gene. Interestingly, in spite of 
amino acid mutation, both S and N protein retained almost 
consistent amino acid properties (especially pH value) 
(Additional file 2: Table $2). Future study will address im- 
portant roles of these polymorphisms in viral replication 
and pathogenesis. In addition, they were divided into two 
distinct small branches (Figs. 1 and 2). These findings sug- 
gested that PDCoVs in southern China have diverged from 
a common ancestor. Despite the emerging genetic diversity, 
overall, PDCoV prevalence is still largely restricted by the 
territory as demonstrated in Figs. 1 and 2. 

For the two enteric coronaviruses (PEDV and TGEV) in 
pigs, the recombination events were often detected. How- 
ever, most of them were from intra-recombination [27-30]. 
Recently, only one emerging recombinant/chimeric virus 
(named swine enteric coronavirus, SeCoV) was discovered 
in swine feces and resulted from inter-recombination of 
PEDV and TGEV, which had a TGEV backbone and a 
PEDV spike gene [31, 32]. In this study, there were no any 
possible recombinant events occurring in PDCoV strains. 
Maybe, the number and length of our obtained PDCoV se- 
quences were limited. In the following study, the recombin- 
ation event of PDCoV warrants further attentions. 


Conclusion 

This study reported the prevalence of PDCoV on swine 
farms in southern China. Phylogenetic analysis of cur- 
rently circulating PDCoV strains in this region and other 
previously reported strains supported the theory of geo- 
graphical clustering of PDCoV infection landscape. The 
origin of various PDCoVs in different countries and re- 
gions should be further studied. 
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